Electronic Structure and Fermiology of PuCoGas 
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By using a relativistic linear augmented-plane-wave method, we clarify energy band structures 
and Fermi surfaces of recently discovered plutonium-based superconductor PuCoGas. We find 
several cylindrical sheets of Fermi surfaces with large volume, very similar to CeMIus (M=Ir and 
Co) isostructural with PuCoGas, in spite of diff'erent /-electron numbers between Ce^^ and Pu^^ 
ions. The similarity is understood by a concept of electron-hole conversion in a tight binding model 
constructed based on the j-j coupling scheme. Based on the present results, we provide a possible 
scenario to explain why a transition temperature is so high as 18. 5K in PuCoGas. 

PACS numbers: 74.25.Jb, 74.70.Tx, 71.18.+y, 71.15.Rf 



Recently it has been discovered that PuCoGas exhibits 
superconductivity [Q. Surprisingly its superconducting 
transition temperature Tc is 18. 5K, which is the highest 
among those yet observed /-electron materials and high 
enough even compared with other well-known intermetal- 
lic compounds. It has been also found that PuRhGas 
becomes superconducting with Tc=8.6K These plu- 
tonium intermetallic compounds PuMGas have the same 
HoCoGa5-type tetragonal structure as CeMIns, a fam- 
ily of cerium-based heavy fermion superconductors 
Note, however, that superconductivity occurs for M=Ir 
(rc=0.4K) and Co (2.3K) in CeMln^, while antiferromag- 
netic (AFM) phase has been found for M==Rh at ambi- 
ent pressure. Another isostructural material including 
uranium is UMGas Q, but superconductivity has not 
been found yet. These HoCoGas-type compounds are 
frequently referred to as "115". 

Regarding superconducting mechanism in the 115 com- 
pounds, first let us consider Ce-115. It has been widely 
considered that it is unconventional d-wave superconduc- 
tor induced by AFM spin fluctuations. In fact, there are 
some evidences such as behavior in nuclear relaxation 
rate p| and node structure measured by thermal conduc- 
tivity ^. For the phase diagram of Ce(Co,Rh,Ir)In5 0, 
AFM phase is found to exist in adjacent to the super- 
conducting phase. Those experimental facts remind us 
of high-Tc cuprates, but a clear difference from cuprates 
should be remarked. Namely, high-Tc superconductivity 
in cuprates occurs by hole doping into AFM insulators, 
while in Ce-115, no hole doping is needed. To under- 
stand the appearance of superconductivity induced by 
AFM spin fluctuations without hole doping, a crucial 
role of orbital degree of freedom has been pointed out 
by Takimoto et al. 

Concerning Pu-115, it is still premature to draw a 
definitive conclusion about the mechanism of supercon- 
ductivity, but we notice that some normal-state proper- 
ties in PuMGas seem to be dominated by AFM spin fluc- 
tuations, e.g., the Curie- Weiss behavior in magnetic sus- 
ceptibility and electric resistivity in proportion to T^-^^ 
m. Thus, it may be natural to consider that supercon- 



ductivity in Pu-based compounds is also induced by AFM 
fluctuations. However, several problems still exist, even 
if d-wave superconductivity is confirmed in both Ce-115 
and Pu-115 materials. One question is, of course, why Tc 
is so high in Pu-115. As is well known, due to difference 
in spatial extension of wavefunctions, 5/ electrons have 
intermediate nature between localized 4/ and itinerant 
Sd electrons. Namely, energy scale of 5/-electrons should 
be larger than that of 4/-electrons, leading to higher Tc 
in 5/ electron systems if we assume the same electronic 
mechanism for superconductivity. However, this cannot 
be the whole story and the situation is not so simple, 
since Pu'^^ ion includes five /-electrons, in contrast to 
one /-electron in Ce'^+ ion. Furthermore one has to ad- 
dress the question why U-115 does not exhibit supercon- 
ductivity. If we follow the above scenario about energy 
scale, U-115 can be superconducting with relatively high 
Tc, but that is not the case. Thus, it is not sufficient to 
consider Pu-115 and U-115 as simple analogues to Ce- 
115 with large energy scale, based only on difference in 
itinerant nature between 4/- and 5/-electrons. 

In this Letter, in order to clarify those points, we 
calculate energy band structures and Fermi surfaces for 
PuCoGas by applying a relativistic linear augmented- 
plane-wave (RLAPW) method. It is found that several 
sheets with large volume form cylindrical Fermi surfaces, 
quite similar to CeMIns (M=Ir and Co) ||l^. The whole 
energy scale in the band structure of PuCoGas is larger 
than that of CeMIns, as naively expected. On the other 
hand, for UCoGas, we have found only several small 
pocket Fermi surfaces to show semi-metal like behavior 
]ll[ , consistent with the fact that UMGas is not super- 
conducting. The similarity between Ce-115 and Pu-115 is 
understood by a simple tight-binding rnodel constructed 
based on the j-j coupling scheme |l2|. A remarkable 
fact is that Pu-115 can be regarded as a hole version of 
Ce-115. Thus, we can conclude that both Ce-115 and 
Pu-115 have the same electronic origin for superconduc- 
tivity, suggesting that Pu-115 has higher Tc due to the 
combination of the electron-hole conversion concept and 
the energy-scale discussion. 



2 




1.0 



PuCoGas 



rr- 0.5 
0.0 



0.5 




r Mxrz ARZX rma 

EYAASTU WV 



FIG. 1: Energy band structure for PuCoGas obtained by the 
RLAPW method. Ef denotes the position of the Fermi level. 




FIG. 2: Calculated Fermi surfaces of PuCoGas for (a) 15th 
band hole sheets, (b) 16th band hole sheets, (c) 17th band 
electron sheets, and (d) 18th band electron sheets. Colors 
indicate the amount of 5/ angular momentum character on 
each Fermi surface sheet and red-shift indicate the increase of 
the / charcter. The center of the Brillouin zone is set at the 
r point. 



First let us briefly explain the RLAPW method. Read- 
ers interested in the formalism can consult with Refs. 
When we calculate the electronic energy band structure 
of 4/ and 5/ compounds, in general, relativistic effects 
should be included, since electrons near the heavy nu- 
cleus must move with a high speed to keep their station- 
ary motion. In order to take into account major rela- 
tivistic effects such as the relativistic energy shifts, the 
relativistic screening effects, and the spin-orbit interac- 
tion, Loucks derived a relativistic augmented-plane-wave 
method based on the Dirac one-electron wave equation 
. Several problems in his method have been improved 
by Hasegawa and co-workers The local density ap- 
proximation is used for the exchange and correlation po- 
tential and spatial shape of one-electron potential is de- 
termined in the mufhn-tin approximation. Self-consistent 
calculations are performed by using the lattice constants 
determined experimentally |l|]. 

In Fig. 1, we show the energy band structure for 
PuCoGas along the symmetry axes in the Brillouin zone 
in the range from — 0.5Ryd. to l.ORyd., where Ryd. indi- 
cates Rydberg and lRyd.=13.6eV. First note that in the 
vicinity of the Fermi level Ep located at 0.446 Ryd., there 
occurs hybridization between Pu 5/ and Ga Ap states. 
Above Ep near the M point, the flat 5/ bands split into 
two groups, corresponding to the total angular momen- 
tum j=5/2 (lower bands) and 7/2 (upper bands). The 
magnitude of the splitting A between the two groups is 
estimated as A(Pu)=l eV, which is almost equal to the 
spin-orbit splitting in the atomic 5/ state of Pu. Note 
that each Pu APW sphere contains about 5.2 electrons 



in the / state, suggesting that valence of plutonium ion 
is Pu'^^, consistent with experimental result Q. 

By using the total density of states at Ep, evaluated 
as N{Ep)—97.3 states /Ryd. cell, the theoretical specific 
heat coefficient 7band is estimated as 16.9 mJ/K^-mol, 
while the experimental electronic specific heat coefficient 
7cxp is 77 mJ/K^-mol [Q. If we define the enhance- 
ment factor for the electronic specific heat coefficient as 
A=7oxp/7band~li wc obtain A=3.6, which is smaller than 



A=10 for CeCoIns |10 . Note that the enhancement of A 
from unity is a measure of electron correlation effect. The 
moderate A in Pu-115 suggests that the correlation effect 
in Pu-115 should be weak compared with Ce-115. Since 
localized nature is stronger in 4/ electrons, the correla- 
tion effect is more significant in Ce-115. 

Now we discuss the Fermi surfaces of PuCoGas. In 
Fig. 1, the lowest fourteen bands are fully occupied. 
The next four bands are partially occupied, while higher 
bands are empty. Then, as shown in Figs. 2(a)-(d), 15th, 
16th, 17th, and 18th bands crossing the Fermi level con- 
struct the hole or electron sheets of the Fermi surfaces, 
summarized as follows: (a) The Fermi surface from the 
15th band includes one small hole sheet centered at the 
r point, (b) The 16th band constructs a large cylindrical 
hole sheet centered at the P point, while two equivalent 
small hole sheets are centered at X points. Depending 
on the energy resolution in the calculation, it is subtle 
whether those two hole sheets touch each other or not, 
but in any case, the quasi two-dimensional large hole 
sheet centered at the P point gives the main contribu- 
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FIG. 3: Energy band structure around Ef for (a) CeCoIns, 
(b) UCoGas, and (c) PuCoGas. In each panel, red and blue 
curves indicate the upper and lower bands to construct the 
hole and electron Fermi-surface sheets, respectively. Namely, 
those are 14th (red) and 15th (blue) for CeCoIns, 15th (red) 
and 16th (blue) for UCoGas, and 16th (red) and 17th (blue) 
for PuCoGas. (d) Calculated Fermi surfaces of CeCoIns. 



tion. (c) The 17th band has a large cylindrical electron 
sheet centered at the M point, (d) The 18th band pro- 
vides another cyhndrical electron sheet centered at the 
M point. 

Let us consider the main Fermi surfaces from the view- 
points of the Fermi-surface volume and /-electron ad- 
mixture. As for the Fermi surface constructed from 18th 
band, /-electrons are not uniformly distributed on it, as 
expressed in color scale. Around the A-point, /-electron 
admixture is large, while p-electron gives a large con- 
tribuation around the M-point. If we ignore three di- 
mensionality and small- volume Fermi surfaces, the main 
contributors are the hole sheet from the 16th band cen- 
tered at the r point and the electron sheet from the 17th 
band centered at the M point. 

In order to gain deeper insight into electronic structure 
of Pu-115, it is instructive to compare with the results 
for Ce-115 @ and U-115 |ll|. In Figs. 3 (a)-(c), we 
show the energy band structures around Ep for CeCoIns, 
UCoGas, and PuCoGas, respectively, in the same en- 
ergy scale. Note that CeCoIns becomes superconducting 
with rc=2.3K j^, while UCoGas is PauH paramagnet [Q. 
First, we can estimate A(Ce)=0.4eV and A(U)=0.8eV, 
which are almost equal to the spin-orbit splittings in the 
atomic 4/ and 5/ states for Ce and U, respectively. As 



expected, we obtain A(Ce)<A(U)<A(Pu). As shown 
in the figure caption, the number to label red and blue 
curves increase one by one in the order of CeCoIns, 
UCoGas, and PuCoGas, corresponding to the increase in 
/-electron number by two per site. Note that shapes of 
red and blue curves among three 115 compounds are sim- 
ilar to one another, since overall band structure around 
the Fermi level is always determined by hybridization 
between broad p-bands and narrow /-bands for 115 com- 
pounds. The center of gravity of the j=5/2 states in 
CeCoIns is about 0.4eV above Ep, while the center of 
those in PuCoGas is slightly lower than Ep ■ Concerning 
UCoGas, the j—5/2 states seem to be just at the Fermi 
energy. This trend is consistent with the number of /- 
electrons in each compounds. The width of j=7/2 and 
5/2 bands around at the M- or A-points becomes broad in 
the order of CeCoIns, UCoGas, and PuCoGas, consistent 
with the difference in 4/- and 5/-electron wavefunctions. 

Here we emphasize that Ce-115 and Pu-115 exhibit 
large Fermi suafces, as shown in Figs. 2 and 4(d). In 
particular, we see a clear similarity between main Fermi 
surfaces of CeCoIns and PuCoGas, except for fine struc- 
tures. Considering only /-electron dominant Fermi sur- 
face with large volume, we observe in common the large 
hole sheet centered at the F point and the large cylin- 
drical electron sheet centered at the M point. On the 
other hand, U-115 has small-pocket Fermi surfaces, as 
deduced from Fig. 2(b) |jri|. Namely, U-115 is consid- 
ered as a semi-metal, which seems to be closely related 
to the reason why U-115 does not exhibit superconduc- 
tivity. The origin of the semi-metallic behavior may be 
traced back to slight overlap among the j—5/2 /-bands 
strongly hybridized with the p-states from Ga ions. 

In order to understand the similarity in energy band 
structures and Fermi surfaces between Ce-115 and Pu- 
115, it is convenient to reanalyze the tight-binding model 
obtained based on the j-j coupling scheme . To con- 
sider the 115 systems, we include only /- and p-electrons 
in the two-dimensional network composed of Ce and In 
(Pu and Ga) ions |l^. Due to the lack of space, we skip 
the details of the model for j=5/2 sextet |l^, but 
the Hamiltonian H is written as H=H{+Hp+Hfp, where 
Hf, Hp, and Hfp are, respectively, /-electron hopping, p- 
electron hopping, and f-p hybridization terms, which are 
characterized by the Slater integrals {ffcr), (ppcr), and 
{fpcr), respectively. Note that crystalline electric field 
terms are simply ignored, since those are much smaller 
than the energies considered here. 

In Fig. 4(a), we show the direct comparison between 
the RLAPW and tight-binding results for (//cr)=4500K, 
(/p(j)=6000K, and {ppa) = 18630K. Top and bottom of 
the tight-binding bands are determined by comparison 
with the RLAPW ones with significant amount of Ga 
Ap states. The Fermi level for H is determined so as 
to include five / electrons. First, overall features of the 
bands in the vicinity of Ep are well reproduced by the 
mixture of broad p- and narrow /-bands. Second, mag- 
nitude of parameters for PuCoGas are large compared 
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FIG. 4: (a)Energy band structures for PuCoGas around Ef 
for the tight-binding model (solid curves) and the RLAPW 
results(dashed curves). Fermi-surface lines discussed here 
are (b)16th band hole sheets, (c)17th band electron sheets, 
and (d)18th band electron sheets. Note that solid and bro- 
ken curves denote the tight-binding and RLAPW Fermi sur- 
faces, respectively, (e) Configurations for /-electrons accom- 
modated in three Kramers doublets. Up and down arrows 
denote pseudo-spin up and down, respectively. 

with those for CeCoIns, (//cr)=4400K, (/pcr)=5360K, 
and (pp(t)=5730K Note that the difference in (ppa) 
between PuCoGas and CeCoIn5 is mainly due to the dif- 
ference of Ga 4p and In 5p electronic states. Then, we 
conclude that 5/ electrons are more itinerant than 4/ 
ones from the present results for tight-binding fitting. 
As shown in Figs. 4(b)-(d), the main Fermi surfaces are 
well reproduced by the tight-binding model. Good agree- 
ments between RLAPW and tight-binding results indi- 
cate the validity of the j-j coupling scheme for PuCoGas. 



In the tight-binding model constructed based on the 
j-j coupling scheme, we are allowed to consider the /" 
configuration with n>2 by accommodating n electrons 
among three Kramers doublets. In Fig. 4(e), we show 
several /-electron configurations. First of all, we note 
the electron-hole conversion relation between Ce"^+ and 
Pu^^^ ions |T^. Thus, Pu-115 can be considered as a 
hole version of Ce-115 and this is the very reason why 
common Fermi surfaces are observed. 

Note that UMGas has been found to be AFM metal for 
M=Ni, Pd, and Pt, while Pauli paramagnetic for M=Fe, 
Co, and Rh Thus, the present band-structure calcula- 
tion assuming the paramagnetic phase is consistent with 
the experimental result for UCoGas. On the other hand, 
a hint to understand AFM metallic behavior for M=Ni, 
Pd, and Pt may be found in the local spin structure. Al- 
though it is difficult to determine the exact valence of 
uranium ion, it should be between U^"*" and U'^^. As 
shown in Fig. 4(e), for U'*+ (U"^^) ion, local spin S—1 
(3/2) may be formed due to the Hund's rule couphng 
and thus, the AFM phase will be favored. However, it is 
still an open problem to explain the metallic behavior in 
the AFM phase as well as the difference in the AFM spin 
structure between UPtGas and UNiGas Q]. 

Finally, we provide one short comment on Np-115. For 
Np'^^ ion, as shown in Fig. 4(e), we can regard it as a hole 
version of U^^ and thus, Np-115 may not exhibit super- 
conductivity, but antiferromagnetism or paramagnetism, 
as an analogue of U- 115. 

In summary, we have performed the band-structure 
calculation for PuCoGas and obtained the Fermi surfaces 
similar to Ce-115 materials. This similarity can be un- 
derstood by the electron-hole conversion picture based 
on the j-j coupling scheme. We believe that high Tc in 
PuCoGas can be understood by combining our electron- 
hole picture with the energy-scale difference in 4/- and 
5/-electrons. 

We thank Y. Onuki, J. L. Sarrao, T. Takimoto, F. 
Wastin, and H. Yamagami for discussions. T. H. and K. 
U. are supported by the Grant-in- Aid for Scientific Re- 
search from Japan Society for the Promotion of Science. 



[1] 
[2] 



[4 



J. L. Sarrao et ai, Nature (London) 420, 297 (2002). 
F. Wastin et al, preprint. 
[3] H. Hegger et al, Phys. Rev. Lett. 84, 4986 (2000); C. 
Petrovic et al, Europhys. Lett. 53, 354 (2001); C. Petro- 
vic et al., J. Phys.: Condens. Matter. 13, L337 (2001). 
Y. N. Grin et al, J. Less. Common Met. 121, 497 (1986); 
Y. Onuki et al, Acta Phys. Pol. B32, 3273 (2001); Y. 
Tokiwa et al, J. Phys. Soc. Jpn. 70, 2982 (2001); Y. 
Tokiwa et al., J. Phys. Soc. Jpn. 71, 845 (2002); H. Kato 
et al, J. Phys. Chem. Solids 63, 1197 (2002). 
[5] Y. Kohori et al, Eur. Phys. J. B18, 601 (2000); G.-q. 

Zheng et al, Phys. Rev. Lett. 86, 4664 (2001). 
[6] K. Izawa et al., Phys. Rev. Lett. 87, 057002 (2001). 
[7] P. G. Pagliuso et al, Physica B312-313, 129 (2002). 



[8] 

[9] 
[10] 

[11] 
[12] 
[13] 



[14] 



T. Takimoto et al, J. Phys.: Condens. Matter. 14, L369 
(2002). 

T. Moriya and K. Ueda, Adv. Phys. 49, 555 (2000). 
T. Maehira et al., J. Ph ys. Soc. Jpn. Supp l. 71, 285 
(2002); T. Maehir a et aL, |cond-mat/021136oi 



T. Maehira et al, |cond-mat/020757E 



T. Hotta and K. Ueda, |cond-mat/02il46'i . 
A. Hasegawa and H. Yamagami, Prog. Theor. Phys. 
Suppl. 108, 27 (1992); H. Yamagami and A. Hasegawa, 
J. Phys. Soc. Jpn. 59, 2426 (1990); M. Higuchi and A. 
Hasegawa, J. Phys. Soc. Jpn. 64, 830 (1995). 
T. L. Loucks, Augmented Plane Wave Method (Benjamin, 
New York, 1967). 



